***
        PROGRAM popbin
    
        implicit none
        integer i,j 
        real *8 f , dum, beta,qnow,ufun
        real :: aaa(10)
        DATA aaa /1.1, 1.5, 2, 3, 4 ,6, 8, 10,14,20/  
        real bbb(10)

        do i=1,11
            do j=1,10
                f=105.0-5.0*i-0.1
                ! print*,(aaa(j))
                beta=3.0
                qnow=15.0
                dum=aaa(j)
                dum=Ufun(f,dum,beta,qnow)
                ! print*,f,aaa(j),dum
                bbb(j)=dum
                
            enddo
            print*,bbb
        enddo

*
      STOP
      END
***
